Tunneling through ferromagnetic barriers on the surface of a topological insulator 
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We study the transmission through single and double ferromagnetic barriers on the surface of 
a topological insulator. By adjusting the gate voltage and magnetization oreintation, the ferro- 
magnetic barrier can be tuned into various transmission regions, where the wavevector-dependent 
tunnelings are quite different. We find that the Klein tunneling can be manipulated or even be 
turned off. These special properties offer the possibility to control electron beams on the "topo- 
logical metal". Various novel devices, such as electronic collimation, wavevector filter, magnetic 
and electric switchs, and wavevector-based spin valve, may be constructed based on our observed 
phenomena. 



I. INTRODUCTION 



Due to its unique electronic properties, the study 
about the graphene system has recently drawn great 
attentions'^. The charge carriers in single-layer 
graphene behave like "relativistic" chiral massless par- 
ticles and possess a gapless linear energy spectrum. Ef- 
fectively, in low energy regime, graphene can be viewed 
as a two-dimensional (2D) Dirac fermion system. For 
such Dirac-like quasiparticles, the transport properties 
are dramatically different from the conventional semi- 
conductors. There is hope to utilize these peculiar trans- 
port characteristics to develop novel relativistic electronic 
devices 4 . 

Recently, another 2D Dirac fermion system has been 
confirmed experimentally. It is called topological 
metal (TM), i.e., the surface states of a strong topolog- 
ical insulator (STI)2^±iJ^±i. The strong topological 
insulator is a new state of matter protected by the time- 
reversal invar iance. Though still with a bulk gap, it is dis- 
tinguished from the normal insulator by the topological 
Z2 invar iance. This distinction results in the existence of 
the gapless electronic states on the sample surfaces. The 
surface state is the TM mentioned above, which can be 
effectively viewed as a 2D chiral Dirac fermion system. 
Due to the time-reversal symmetry and the topology of 
the bulk gap, the topological surface states are robust 
against disorder. 

Because of the similar energy dispersions, the trans- 
port property of the TM is very similar to that of a 
graphene. However, there are important differences be- 
tween the two systems. TM is a surface state of semicon- 
ductor alloy, e. g., Bii- x Sb x or Bi2Se^ 1 thus more stable 
than a fragile thin graphene sheet. Furthermore, the spin 
of the Dirac fermion of a TM is the real spin, while for 
graphene it is the pseudo spin. Hence, the magnetic or 
spin control of the electronic transport on a TM can be 
much more effective than that on a graphene. In this 
work, we study the transport properties with ferromag- 
netic barriers on top of a topological surface. As shown 




FIG. 1: Schematics of the ferromagnetic barriers on topolog- 
ical surface. Double barriers are shown as an example. The 
whole system has been divided into five regions. The electro- 
static potential of region i is V% (i = 1 ~ 5). U denote the 
interface between region i and i + 1. 



in Fig.l, for a ferromagnetic barrier, a ferromagnetic in- 
sulator is put on the top of the topological insulator to 
induce an exchange field via the magnetic proximity ef- 
fect. Meanwhile, a metal gate on the top of the insulator 
is used to control the electrostatic potential. We find 
that by adjusting the magnetization orientation in the 
ferromagnetic insulator and the gate voltage, the ferro- 
magnetic barrier can be tuned into various transmission 
regions. The behaviors of the transmission probability, 
as a function of the incident angle of the electron, in 
different regions are rather different. An interesting is- 
sue is about the Klein tunneling. It is well known that, 
for Dirac fermion, merely with electrostatic potential one 
cannot block the electron transport, and the normal in- 
cidence is at least partially transmitted. But for the fer- 
romagnetic barrier, because the induced exchange field, 
not only could the angle of Klein tunneling be controlled, 
but also the Klein tunneling could even be blocked. We 
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carefully analyze the transmission probabilities in vari- 
ous regions, for both single and double barriers. Based 
on these analyses, we propose that electron beam collima- 
tion, wavevector filtering, magnetic and electric switch, 
wavevector-based in-plane spin valve could be realized in 
such structures. These transport phenomena are quite 
different from the case with ferromagnetic barrier on 
graphen o 15 i 16 i 17 i 18 . 

This paper is organized as follows. In Sec. II, we in- 
troduce the model and formalisms. In Sec. Ill, we dis- 
cuss the analytical and numerical results of the transport 
properties of single and double ferromagnetic barriers. 
And we then show how to construct novel electronic de- 
vices based on the phenomena obtained in our theory. 
Finally, a brief summary is given in Sec. IV. 



II. MODEL AND FORMALISM 

The Dirac fermion of the topological surface can be 
described as: 



H = hv f k- a + M-(T 



(1) 



where M = M(sin rn cos </> m , sin rn sin </> m , cos rn ) is the 
effective exchange field induced via magnetic proximity 
effect. In this work, we only consider the in-plane ex- 
change field. The main consequence of the z compo- 
nent of the exchange field is to open an energy gap and 
its effects have been investigated in a recent theoretical 
study 19 . If the Fermi level of the whole system is in- 
side the gap, the transport of the barrier is completely 
blocked. No tunneling is permitted. Otherwise, the ef- 
fect is just a change of the position of the Fermi surface. 
But in our ferromagnetic barrier model, adjusting the 
gate voltage will also change the Fermi surface. So for 
simplicity and without loss of generality, we merely con- 
sider an in-plane exchange field. The sketch of the barrier 
structure is shown in Fig.l. 

Take the single barrier structure for example, the whole 
system consists of three different regions and the inter- 
faces between these regions are labeled by k. With the 
above Hamiltonian, the eigenfunctions for region a are 
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wave functions at the interface l\ = 0, we get the transfer 
matrix 
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The transfer matrix at the interface I2 = a can be repre- 
sented as 
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where is the same as except the region index. 
Finally, we will get the total transfer matrix for the single 
barrier 



T s = T h T h 



(6) 



With the transfer matrix, the reflection and translation 
coefficients can be got through following equation 



w; 



(7) 



T = \t\ 2 is just the transmission probability we want. 

Double barrier structure could be treated in similar 
way. The only difference is that we have four interfaces 
here. 



rri rri rri rri rri 

J-d — 1 lA 1 h 1 l2 1 li 



(8) 



With the transmission probability, in the linear trans- 
port regime and for low temperature, we can get the con- 
ductance density (the conductance per unit length along 
the interface) of the barrier structure^ 



,7T/2 

g = g dO in T(e f , in ) cos Q { , 

J-7T/2 



(9) 



Here, go = 4 ^2^i Vf , is the Fermi energy of the system, 
9 i n is the incident angle of the electron and T(e/,^ n ) is 
the transmission probability through the barrier. 



where S a = e - V a: e z9a = (k ax + ik ay )/\k a \, and k^ x = 

-Mcos(/w ± V( e ~ V a) 2 ~ V^y + Msin0 ma ) 2 . V a is 
the potential of this region, e is the energy of the incident 
electron. Due to the translation invariant along the y- 
axis, the momentum k y is conserved, i.e. k ay = k y . 

In order to solve the transport problem, we use the 
standard transfer matrix method^. By matching the 



III. RESULT AND DISCUSSION 

A. Single barrier 

For a single barrier, an analytical expression of the 
transmission probability could be given. 
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FIG. 2: (Color online) The overlap of the Fermi surfaces of 
region 1 and 2. Here, region 2 is a ferromagnetic barrier 
and its fermi surface (red circle) has been shifted due to the 
induced exchange field M. e is the energy of incident electron. 
Vi(2) is the electrostatic potential of region 1(2). m is the 
in-plane rotational angle of the exchange field. 



where q x = k^ + M cos (j) m . Formally, it is the same as 
the one for the graphene potential barrier except a shift 
of the Fermi surface due to the in-plane exchange field. 

As mentioned above, for a ferromagnetic barrier, 
changing the gate voltage will alter the Fermi surface, the 
in-plane exchange field will also cause a shift. The key 
issue about the transport is that, during the tunneling, 
the energy and k y of the incident electron are conserved. 
As shown in Fig. 2, in the incident region (Region 1 in 
Fig.l), the range of k y values |e — V\\ > \k y \ is plotted. 
But not all the incident electrons are allowed to trans- 
mit, since the Fermi surface has been changed in the fer- 
romagnetic barrier ( Region 2 in Fig. 1). The k y range 
of the ferromagnetic barrier is |e — V2 1 > | k y + M sin m | . 
Therefore, the permitted incident angle of the electron is 
determined by the overlap of the k y ranges of adjacent re- 
gions. It means that, through adjusting the gate voltage 
and the direction of the magnetization of the ferromag- 
netic barrier, we can tune the acceptable range of k yi and 
rather different transmissions will appear accordingly. 

In the case of Fig. 3, V2 is far from the energy of 
incident electron e. So for any in-plane rotation angle 
(j) m , there are always | e — V2 1 — M sin (j) m > \ e — V\ | and 
— |e — Vi I > — I e — V2 1 — M sin m . Hence, all the incident 
angles are allowed from the viewpoint of k y conservation. 
We can say that the ferromagnetic barrier is in the total 
transmission status now. As shown clearly, when we ro- 
tate the magnetization orientation, the angle of the Klein 
tunneling will be changed and new resonance peaks will 
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FIG. 3: (Color online) The transmission probability of the 
single ferromagnetic barrier in the total transmission status. 
Bin is the incident angle of the electron. m is the in-plane 
rotation angle of the magnetization of the barrier. Here, Vf — 
10 6 ras _1 , e 80meV, Vi = 0, V 2 = 200raeF, M 40meV, 
h =0A 2 = lOOnra. 
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FIG. 4: (Color online) The transmission probability of the 
single ferromagnetic barrier in the partial transmission status. 
V2 = 90meV and other parameters are the same as Fig. 3. 



appear. 

Another intriguing case occurs when the gate volt- 
age of the barrier V2 is near the energy of the incident 
electron e. In this situation, the relations, |e — V\\ > 
\e-V 2 \ -Msm(l) m and -\e-V 2 \ -Msin0 m > -|e- Vi\ 
are always hold for any angle </> m . As shown in Fig. 
4, we can get a transmission sector of incident angle in 
this partial transmission status. The width of the sector 
is determined by |e — V2I and its position can be tuned 
through adjusting the magnetization orientation. It obvi- 
ously can be used as a wavevector filter. Furthermore, as 
pointed in electron optics, this has immediate applicabil- 
ity in electron beam collimation and lensin g 21 i 22 . Better 
collimation effect could be achieved through a series of 
barriers. 

When —M — |e — V2 1 > |e — Vi\, for some angle m , 
the overlap between the k y ranges of the two regions will 
be absent and the conductance becomes zero. This is 



FIG. 5: (Color online) Up:The transmission probability of 
the single ferromagnetic barrier in the blockable status. V2 = 
90meV , V± = 60meV and other parameters are the same as 
Fig. 3. Down: The corresponding conductance density in this 
case. 




FIG. 6: (Color online) Up:The transmission probability of the 
single ferromagnetic barrier in the blockable status. m = 
0.57T and other parameters are the same as Fig. 5. Down: 
The corresponding conductance density in this case. 



the blockade status, in which the transmission of the 
ferromagnetic barrier can be completely blocked. Since 
that both V2 and (j) m could be used to tune the over- 
lap, we can construct magnetic and electric on-off with 
the single ferromagnetic barrier structure. Fig. 5 gives 
the case of magnetic switch. With a definite V2, there 
is a blocked sector of angle m , in which the transmis- 
sion probability is zero for any incident electron. The 



FIG. 7: (Color online) Up:The transmission probability of the 
double ferromagnetic barrier structure. (j)mi{2) is the in-plane 
rotation angle of ferromagnetic barrier 1(2), and (j) m 2 — 0.57T. 
Both barriers are in the partial transmission status. Vi = 
V 3 = V 5 = 60meV, V 2 = V 4 = 90meV, e SOmeV and 
h — h = h — h = U — h = lOOnra. Down: The corresponding 
conductance density in this case 



conductance of the magnetic switch is shown as well. 
The alteration of electrostatic potential can also com- 
pletely block the transmission. This is quite different 
from a normal tunneling barrier for Dirac fermions. It 
is known, due to the Klein paradox, merely with elec- 
trostatic potential one can not block the transport of 2D 
Dirac fermions. However, with a ferromagnetic barrier 
the condition is different. For a definite angle m , if 
|Msin0 m | > |e — Vi\, there must be a critical gate volt- 
age V2C that — Msin^ m — \e — V2 C \ = |e — Vi\. Compared 
with V20 when |e — V2I decreases, the k y ranges of the 
two regions separate out and the transport is turned off 
(See in Fig. 6). It is just an electric switch which the 
on/off state of the current is controlled by the gate volt- 
age. Interestingly, when increasing V2, the current will 
be turned on (See in Fig. 6). This unique property, i.e. 
the voltage window for the open status, is not possible 
with a normal semiconductor barrier. Clearly, this trans- 
port phenomenon could be used to distinguish TM from 
the normal semiconductor and the graphene systems. It 
means that transport measurement is also a very useful 
tool to experimentally confirm the existence of the TM. 



B. Double barrier 

For the double barrier structure, We also could under- 
stand the transport properties through the conservation 
of the k y and the energy. As mentioned above, depending 
on the voltage and magnetization, single ferromagnetic 
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FIG. 8: (Color online) The transmission probability of the 
double ferromagnetic barrier structure. Both barriers are in 
the total transmission status. Vi = V3 = V5 = 0, V2 = V4 = 
200meV, e 80meV and h - h = h - h = U - h = lOOnra. 

barrier could be tuned into three different transmission 
statuses: total transmission , partial transmission , and 
the blockade situations. The transport behavior of the 
double barrier system is more complex than that of the 
single one, since that the two barriers can be tuned into 
different transport situations. 

Here, we investigate two interesting examples. One is 
the case that two ferromagnetic barriers are both in the 
partial transmission status. In this case, both barriers 
have their own transmission sector of the incident angle 
of the electron. Only when the two transmission sectors 
overlap, the double barrier structure becomes conduct- 
ing. Hence, as shown in Fig. 7, the direction of magneti- 
zation of the second barrier is fixed <j) m 2 = tt/2. Only if 
the angle of the first barrier m i is tuned in parallel, the 
transmission of the double barrier structure is nonzero. 
The corresponding conductance is also given in Fig. 7. 
Essentially, in this situation the double barrier system be- 
comes a special wavevector-dependent in-plane spin valve 
system, for which the parallel configuration is the only 
conducting state. The blockade of the transmission is 
mainly because the mismatch of the wavevectors, while 
for normal spin valve the primary reason is the different 
spin dependent potentials. Therefore, for a normal spin 
valve, the conductance will vary continuously with the 
angle between the magnetizations. But, as shown in Fig. 
7, for the in-plane spin valve on TM the change of the 
conductance is rather sharp and the the conductance of 
the non-parallel configuration is exactly zero. These spe- 
cial transport properties of the double barrier structure 
on TM could also be tested in transport experiments as a 
distinction between the TM and normal semiconductors. 

For a single barrier, it is hard to directly observe 
the partial transmission region through the conductance. 
This is because that the change of position of the trans- 



mission sector will not induce an obvious change in the 
conductance. But in the double barrier system this can 
be done. 

As a comparison (See in Fig. 8), we also calculate the 
case that both barriers are in the total transmission sta- 
tus. In this situation, the angle dependent transmission 
probability is much more complicated. 

IV. SUMMARY 

In conclusion, we have investigate the transport prop- 
erties of the ferromagnetic barriers on the surface of 
a topological insulator. The single and double barrier 
structures have been carefully considered. We show that 
the gate voltage and the magnetization could be used 
to manipulate the Fermi surface of the barrier. Due to 
the conservation of energy and k y of an incident electron 
during the transport process, adjusting the Fermi surface 
will tune the single ferromagnetic barrier into three dif- 
ferent transmission statuses: total transmission , partial 
transmission , and the transmission blockade. In the to- 
tal transmission case, all the incident angles are allowed. 
The angle of the Klein tunneling and position of the 
resonance vary with the direction of the magnetization, 
and meanwhile new resonance peaks appear. In the par- 
tial transmission situation, only the incident angles in a 
transmission sector are allowed. The position and width 
of the transmission sector can be adjusted through the 
gate voltage and the magnetization. Therefore, the elec- 
tron collimation and wavevector filter could be realized 
based on this phenomenon. With proper gate voltage and 
magnetization, one can even separate the Fermi surfaces, 
reaching the blockade status. In this case, no transport 
can occur, even the Klein tunneling. We propose that 
both magnetic and electric switchs, based on the single 
ferromagnetic barrier on topological surface, could be re- 
alized in this situation. We also study the transport of 
the double barrier structure. In this system, we can con- 
struct a special wavevector-dependent spin valve when 
both barriers are in the partial transmission status. The 
spin valve is conducting only if the magnetization of the 
two barriers are parallel. We emphasis that these unique 
transport properties of TM could be used to distinguish 
TM from the normal semiconductors and the graphene 
systems. Hence, the observations of these transport prop- 
erties provide confirmation of the existence of TM. 
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